Quasinormal modes and stability of the rotating acoustic black hole: numerical 

analysis 



Vitor Cardoscfl 

Centra de Fisica Computacional, Universidade de Coimbra, P-3004-516 Coimbra, Portugal 

Jose P. S. Lemo^j] 

Centro Multidisciplinar de Astrofisica - CENTRA, 
Departamento de Fisica, Institute/ Superior Tecnico, 
Universidade Tecnica de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

Shijun Yoshid£0 

Science and Engineering, Waseda University, Okubo, Shinjuku, Tokyo 169-8555, Japan 

(Dated: February 7, 2008) 

The study of the quasinormal modes (QNMs) of the 2+1 dimensional rotating draining bathtub 
acoustic black hole, the closest analogue found so far to the Kerr black hole, is performed. Both 
the real and imaginary parts of the quasinormal (QN) frequencies as a function of the rotation 
parameter B are found through a full non-linear numerical analysis. Since there is no change in 
sign in the imaginary part of the frequency as B is increased we conclude that the 2+1 dimensional 
rotating draining bathtub acoustic black hole is stable against small perturbations. 

PACS numbers: 04.70.-s 



I. INTRODUCTION 

A guitar string in vacuum vibrates in normal modes. 
However, a guitar string immersed in air vibrates truly 
in quasinormal modes (QNMs) since it looses energy to 
sound waves. The same happens to an acoustic black 
hole. Due to energy losses to the medium in sound wave 
acoustic perturbations, its modes of vibration are quasi- 
normal (QN). This QN behavior is expected, as acoustic 
black holes are real analogues of true general relativistic 
black holes which are known to vibrate in a QN form due 
to losses in gravitational and other types of radiation to 
infinity. 

Acoustic black holes and their analogies with true 
black holes were found by Unruh by noticing that the 
equations of motion for sound waves in a background 
fluid flow and the equations of motions for a scalar field 
in a black hole background are similar [lj. The back- 
ground fluid flow acts as an effective black hole metric 
implying the existence of a sound horizon, the surface 
where the velocity of the flow is equal to the velocity of 
sound in the medium. A number of phenomena that oc- 
curs with true black holes also occurs with acoustic black 
holes. For instance, Hawking radiation is now phonon 
radiation @, 01 BL Geodesic and causal structure can 
also be studied as was done in [|| where Penrose-Carter 
diagrams for several effective acoustic spacetimes were 
drawn. Black holes in other analogue systems, like con- 
densed matter ones, can also be studied (see 0, to 
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name a few). 

In the study of true general relativistic black holes 
QNMs are important for a number of reasons. The most 
important is that they provide means to identifying black 
hole parameters like the mass and angular momentum, 
since the real and imaginary part of the QN frequencies, 
which give the vibration frequency and the exponential 
damping frequency respectively, depend only on those pa- 
rameters. Through their study one further tests the sta- 
bility of the system, as any imaginary frequency with the 
wrong sign would mean an exponentially growing mode, 
rather than damping. QNMs have also been connected to 
the quantization of the black hole area, where it seems 
that the highly damped modes, those that are almost 
instantaneous, are associated with transitions between 
area levels at large quantum numbers (see |lfj| for a re- 
view on QNMs and a list of complete citations therein). 
Since QNMs appear naturally in general relativistic black 
holes, they also should appear in their analogues, the 
acoustic black holes. 

The study of QNMs of the 2+1 dimensional rotating 
draining bathtub acoustic black hole [J, ll2j , the closest 
analogue found so far to the Kerr black hole (see ElJ 
for developments on this subject), was initiated in |l3| 
through the use of a WKB approximation, valid for small 
black hole rotation. Various interesting points besides 
QNM behavior, like late-time tails and superradiance, 
were also discussed. In relation to large damped QNMs 
it was found that there are no asymptotic QN frequen- 
cies, a puzzling result if connected to the area quantiza- 
tion issue. Moreover, the WKB scheme used in [lj| gave 
indications that, due to a possible change in the sign of 
the imaginary part of the frequency u>, this 2+1 acous- 
tic black hole could become unstable at large rotation 
parameter B, or at large azimuthal number m. Recent 
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work on this [1J], making use of a matching procedure, is 
also not appropriate to give a convincing stability proof. 
Is the Imw part lower, equal or greater than zero for B or 
m large? In this paper we make a full numerical analysis 
to settle this issue. We show that the instability never 
sets in. 



II. FORMALISM 

The draining bathtub model is a 2 + 1 dimensional flow 
with a singular vortex-like sink (or source) at the origin. 
The acoustic metric describing the propagation of sound 
waves in this fluid flow is £|: 



ds 2 



c 
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-drdt - 2Bd<pdt + dr 2 + r 2 



(1) 



Here A and B are arbitrary real positive constants related 
to the radial and angular components of the background 
fluid velocity: 
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where r and <fi are orthogonal unit basis vectors point- 
ing along the axes. It is however better to work with 
a more transparent metric. Some physical properties of 
our draining bathtub metric are more apparent if we cast 
the metric in a Kerr-like form performing the following 
coordinate transformation (see [12. H^): 



dt 
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Then the effective metric takes the form 



ds 2 



A 2 + B 2 



1 - 



A 2 



c 2 dP 



dr 2 - 2Bd<j)dt 



(5) 

As explained in |13| , this metric and the Kerr metric dif- 
fer in an important aspect, in that whereas the rotation 
for the Kerr black hole is bounded from above, here it is 
not, at least in principle. Thus, B could be as large as 
desired. 

The propagation of a sound wave in a barotropic in- 
viscid fluid with irrotational flow, which is assumed to 
be the case, is described by the Klein-Gordon equation 
V M V M $ = for a massless field ^ in a Lorentzian acous- 
tic geometry. Separating variables by the substitution 
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implies that vE'(r) obeys the wave equation 

*,f,f, +Q* = 0, 
where the generalized potential is given by 
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The tortoise coordinate is defined by 

dr* 



dr 



where A = (1 — 1/f 



Explicitly, 



1 



r* = f + - log 



f - 1 



f + 1 



(8) 
(9) 

(10) 



We have also performed the following rescaling: f = 
rA/c, Cj — uiA/c 2 , B = B/A. The rescaling effectively 
sets A = c = 1 in the original wave equation, and picks 
units such that the acoustic horizon fu = 1. From now 
on we shall omit hats in all quantities. The rescaled wave 
equation (JJJ will be the starting point of our analysis of 
QNMs. The characteristic QNMs of the rotating acous- 
tic black hole can be defined in the usual way, imposing 
appropriate boundary conditions and solving the corre- 
sponding eigenvalue problem. Close to the sound horizon 
we seek solutions of equation Q behaving as 



— Br. 



(11) 



Classically, only ingoing waves - that is, waves falling into 
the acoustic black hole - should be present at the horizon. 
This means (according to our conventions on the time 
dependence of the perturbations) that we must choose 
the minus sign in the exponential. At spatial infinity the 
solutions of Q behave as 

* - e w * . (12) 

In this case we require that only outgoing waves (waves 
leaving the domain under study) should be present, and 
correspondingly choose the plus sign in the exponential. 
This boundary condition at infinity may be cause for ob- 
jections. Indeed, no actual physical apparatus will be ac- 
curately described by these boundary conditions: a real 
acoustic black hole experiment will certainly not extend 
out to infinity. However, we may imagine using some 
absorbing device to simulate the "purely outgoing" wave 
conditions at infinity (for another example in which an 
absorbing device modelin g sp atial infinity could be re- 
quired, cf. Section XI of [l5| - in particular their Fig. 
5). 

The boundary conditions Ijllll and H12|) are satisfied 
by only a discrete number of frequencies u>, the QN fre- 
quencies wqn. The QN frequencies are in general com- 
plex numbers, and the imaginary part is usually nega- 
tive, which means that perturbations die exponentially 
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as time goes on (recall that the time dependence of the 
field is e~~ lu)t , and a negative imaginary part for ui means 
exponential decay with time). If all the QN frequencies 
have a negative imaginary part, we say that the system is 
stable. Note also that, for a given to there will in general 
be an infinite number of ujqn satisfying the boundary 
conditions. We shall order these ujqn by imaginary part: 
the ujqn having smallest (in magnitude) imaginary part 
will be called the fundamental (n = 0) frequency, the one 
having the second lowest imaginary part will be the first 
overtone n = 1, and so on. 



four-term recurrence relation to the three-term one, given 
by 

a' ai + 0' o a o = , 

a' k a k +i + I3' k a k + 7^-1 = , ft = 1, 2, • • • , (17) 

where a' k , (3' k , and j' k are given in terms of ak, Pk, 7fe and 
4 by 



a' k = a k , P'k=0k, l'k=lk, for ft = 0,1, (18) 



and 



III. NUMERICAL PROCEDURE 

A. Frobenius expansion 

In order to numerically obtain the QN frequencies, 
in the present investigation, we make use of Leaver's 
method 0, EJ]> which is known to yield excellent re- 
sults. 

Let us introduce a new independent radial variable x, 



defined as x — 



In terms of the new variable x. the 



tortoise coordinate is then reduced to 
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for ft > 2 . 



Now that we have the three-term recurrence relation for 
determining the expansion coefficients a k , the conver- 
gence condition for the expansion 114|) , namely the QNM 
conditions, can be written in terms of the continued frac- 
tion as 



a o7i «i7 2 a 273 

ft" 03-' 
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The perturbation function \l/ may be expanded around 
the horizon as 



1 



1 



— i(oJ— Bm)/2 00 

^a fe (l-x) fe , (14) 

k=0 



where a is taken to be ao = 1. The expansion coeffi- 
cients ak in equation l|14|) are determined via the four- 
term recurrence relation (it's just a matter of substituting 
expression (|14f) in the wave equation 0), given by 



a a ai + f3 a = , 
aia 2 + 0iai + 71 a = , 
a k a k+ i + I3 k a k + "fkO-k-i + 5 k a k -2 
where 



(15) 



0, ft = 2,3, 



a k = -8(1 + k) (1 + ft + iBm - iuj) , 

P k = 4{1 + 5fc 2 + m 2 + ft (5 + AiBm - 8iu) 



Ik 



-4iuj - Auj 2 + 2Bm (i + 2uj)} , 
2(3-8fc 2 -Uk (Bm-2uj)), 
-3 - 4ft + 4fc 2 . 



Since the asymptotic form of the perturbations as — > 
00 (r* — > —00) is written in terms of the variable x as 

giwr, ^ e iuJX ± ( e —i(u—Bm)r, ^ M _ x \—i(u—Bm)/2\ 

the expanded perturbation function ^ defined by equa- 
tion l |l 4(1 automatically satisfy the QNM boundary con- 
ditions if the power series converges for < x < 1 . Mak- 
ing use of a Gaussian elimination [16( , we can reduce the 



where the first equality is a notational definition com- 
monly used in the literature for infinite continued frac- 
tions. Here we shall adopt such a convention. 

An analysis of the large ft behavior of the expansion 
coefficients ak shows 



a k +\ 
iim 

k^oo ak 



1± 



k^ 



Aicu 



4ft 



(21) 



The series expansion (1141) will converge uniformly only 
if the sign of the second term in the right hand side of 
((21(1 is chosen such that Re(±(— 2iui)? ) < 0, which will 
happen only for the QN frequencies u>. In case that the 
convergency of the continued fraction ((20(1 is not very 
good, one can use Nollert's technique to avoid this dif- 
ficulty of the convergence In Nollert's technique, 
equation l(2T|) plays an essential role. 



IV. NUMERICAL RESULTS 



Our numerical results, which are all consistent with 
the results in ^3j> are shown in Figs. 11181 

(i) m > 0: In Fig. Ill II we show results pertaining to 
perturbations having positive to, i.e., co-rotating waves. 
In Figs. Ill2l we show the real part of the QN frequencies 
for to = 1 and to = 2 modes respectively, as a func- 
tion of the black hole rotation. Higher m modes follow 
a similar pattern. One can see from these plots that for 
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FIG. 1: The real part of the QN frequency as a function of 
the rotation parameter B/A, for several overtones of a m = 1 
mode. Here, th = A/c is the horizon radius. Note how all 
the several lowest overtones "coalesce" in the high rotation 
regime, all growing linearly with B/A. 



FIG. 3: The imaginary part of the QN frequency as a function 
of the rotation parameter B/A, for several overtones of a m = 
1 mode. It is clear from this plot that the imaginary part of 
the QN frequencies of m > modes is very insensitive to the 
rotation of the black hole. 




100 



B/A B/A 



FIG. 2: The real part of the QN frequency as a function of 
the rotation parameter B/A, for several overtones of a m — 2 
mode. Again, all different overtones have the same behavior 
for high rotation. 



FIG. 4: The imaginary part of the QN frequency as a function 
of the rotation parameter B/A, for several overtones of a m = 
2 mode. 



low black hole rotation parameter B the different over- 
tones are clearly distinguished, but that as the rotation 
increases they tend to cluster and behave very similarly. 
For very large rotation B, all the overtones behave in the 
same manner, and in this high rotation regime the real 
part of the QN frequency scales linearly with the rota- 
tion. Indeed we find that the slope is also proportional 
to to so that 

TTlBc^ 

Re[wQAr] ~ — jj— as B — > oo , for to > (22) 



We notice that this behavior was already present in the 
WKB investigation in 1.3!]. In Figs. 13141 we show the 
imaginary part of the QN frequencies as a function of the 
rotation parameter, for to = 1 and to = 2 respectively. 
Different overtones have different imaginary parts. Note 
also that for high B the real part of the modes coalesce 
whereas the imaginary part does not. The magnitude of 
Im[a;QAr] increases with B, which was observed also in the 
WKB approach Thus, as the rotation increases the 
perturbation dies off quicker. This also means that the 
black hole is stable against to > perturbations, because 
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FIG. 5: The real part of the QN frequency as a function of the 
rotation parameter B/A, for several overtones of a m = — 1 
mode. Notice that for each overtone number n there is a 
critical rotation at which the mode crosses the axis, i.e., there 
is a critical rotation B/A at which the real part of the QN 
frequency is zero. Higher overtones cross the axis at a slower 
rotation. We have not been able to follow the mode beyond 
this point. 



FIG. 7: The imaginary part of the QN frequency as a func- 
tion of the rotation parameter B/A, for several overtones of 
a m — — 1 mode. We have not been able to follow the modes 
beyond a certain critical point (defined as the rotation B/A 
for which the real part of the QN frequency is zero- see figs. 
15161) . Nevertheless, an judging by the modes we did manage to 
follow, namely the fundamental mode, it seems that Im[wQjv] 
never crosses the axis, i.e., it is always negative, and therefore 
the mode is stable. 




FIG. 6: The real part of the QN frequency as a function of the 
rotation parameter B/A, for several overtones of a m — —2 
mode. 



FIG. 8: The imaginary part of the QN frequency as a function 
of the rotation parameter B/A, for several overtones of a m = 
—2 mode. 



the imaginary part is always negative. 

(ii) m < 0: In Figs. 15181 we show results concerning 
perturbations having negative m, i.e., counter-rotating 
waves. The behavior of the QN frequencies for m < 
is drastically different from the m > perturbations. In 
Figs. 15161 we plot the dependence of Re[wQjv] as a func- 
tion of the rotation of the black hole B. As B increases 
the magnitude of the real part of the QN frequency de- 



creases. The oscillation frequencies for the fundamental 
modes, labeled by n = 0, indeed get close to the hori- 
zontal axis as B goes to infinity. However, we haven't 
been able to track some overtone modes with negative m 
for very high rotation since, as can be seen in Figs. 15161 
the real part of these modes eventually change sign. It 
is extremely difficult, using the method employed here, 
to compute modes having RefwQAr] ~ 0. Nevertheless, 
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supposing that (as the numerical studies for the funda- 
mental modes indicate) the QN frequencies asymptote 
to zero for very large B, a WKB [2£j analysis shows that 
ujqn ~ k/B, where k is some m-dependent constant. 
The imaginary part of the QN frequencies behaves in a 
similar manner, as seen in Figs. 1 7181 
(iii) 772 = 0; For circularly symmetric (m = 0) modes, 
our numerical method shows no sign of convergence. For 
777 = 0, the equation (JJJ can be written in the simpler 
form 



where 



V 



f 2 - 1 



1 5 

4-7*2 ^4 



(23) 



(24) 



The potential V is not positive definite, and this pre- 
cludes also a simple stability proof. 



V. CONCLUSIONS 

In this paper we have studied numerically the quasinor- 
mal modes of the (2 + l)-dimensional draining bathtub 
metric, which describes a rotating acoustic black hole. 
Our results indicate that this is a metric stable against 
small perturbations, although one would like to have also 
an analytical proof of this statement. Notice that this 
proof, would most certainly encompass also the the sta- 
bility of the usual general relativistic black holes against 



perturbations of a charged scalar field; In fact, for such 
a field it is possible to show |2l| that the generalized po- 
tential goes like (ui — eq/r) 2 — V, where e is the scalar 
field charge, q the charge of the black hole, and V an 
^-independent potential. This is of the same form as the 
generalized potential dealt with here, equations (|23|) -lj%)l. 
Ever since the first studies on the Kerr geometry (see for 
example p ^ . and for numerical results regarding the 
QN frequencies of rotating Kerr black holes), one knows 
that it is extremely difficult to prove stability when uj 2 — Q 
is not positive definite (or when it is w-dependent), and 
therefore the special case Q = (lo — eq/r) 2 — V could 
shed some light on this problem. To conclude, we would 
like to draw attention to the fact that the results pre- 
sented here are similar in many respects to the results 
concerning higher dimensional rotating black holes |24| . 
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